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ABSTRACT 

In a virialized stellar system, the mean-square velocity is a direct tracer of the energy 
per unit mass of the system. Here, we exploit this to estimate and compare root- 
mean-square velocities for a large sample of nuclear star clusters and their host (late- 
or early-type) galaxies. Traditional observables, such as the radial surface brightness 
and second-order velocity moment profiles, are subject to short-term variations due 
to individual episodes of matter infall and/or star formation. The total mass, energy 
and angular momentum, on the other hand, are approximately conserved. Thus, the 
total energy and angular momentum more directly probe the formation of galaxies 
and their nuclear star clusters, by offering access to more fundamental properties of 
the nuclear cluster-galaxy system than traditional observables. We find that there is 
a strong correlation, in fact a near equality, between the root-mean-square velocity 
of a nuclear star cluster and that of its host. Thus, the energy per unit mass of a 
nuclear star cluster is always comparable to that of its host galaxy. We interpret 
this as evidence that nuclear star clusters do not form independently of their host 
galaxies, but rather that their formation and subsequent evolution are coupled. We 
discuss how our results can potentially be used to offer a clear and observationally 
testable prediction to distinguish between the different nuclear star cluster formation 
scenarios, and/or quantify their relative contributions. 

Key words: galaxies: nuclei - galaxies: formation - galaxies: photometry - methods: 
statistical - galaxies: kinematics and dynamics - galaxies: evolution. 


1 INTRODUCTION 

Recent studies have revealed that nuclear star clus¬ 
ters (NSCs) occu r commonly in galaxies o f ev¬ 
ery type (e.R. iB inggel i, Sandage fc Tarenghil 1 19841 : 
iKormendv Sz Kennicutin2004h These massive and com¬ 
pact star clusters occupying the centres of galaxies are 
characterized by small effective radii and high central 
surface brightnesses that help to distinguish them from the 
surrounding bulge and disk. 

Several theories have been proposed to explain the ori¬ 
gins of NSCs, but two in particular have emerged as leading 
scenarios. The first invoke s gas accretion at the centres of 
galaxies to form stars (e.g. [McLaughlin. King fc NavakshinI 
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C.Knigge@soton.ac.uk (CK); iskren.y.g@gmail.com (IG); den- 
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I2OO6I : iNavakshin. Wilkinson fc King||200^ . The second pos¬ 
tulates successive mergers of globular clusters (GCs) that 
spir al into the galactic centre due to dynamical friction 


(e.g. [Tremaine. Ostriker fc Spitz^ll975l: Quinlan fc Shapir3 


Il99d : Gnedin. Ostriker fc Tremaine! 2014ll . Several authors 
have used (primarily N-body) simulations to succesfully re¬ 
produce some of the observed features of NSGs formed via 
this mechanism, includi ng their effective radii and central ve - 
locity dispers ions (e.g. Gapuzzo-Dolcett a fc Miocchi|[200al . 
For example, iGnedin. Ostrike~fc Tremaine! ( 2014h recently 
modeled the effects of dynamical friction on the orbits of 
GGs within their host galaxies, in conjunction with the un¬ 
derlying dynamical evolution of the infalling clusters, in¬ 
cluding mass loss (due to both stellar evolution and tidal 
stripping) and tidal disruption. The model successfully re¬ 
produces the observed features of NSCs in both late- and 
early-type galaxies. 


Whatever the dominant formation scenario for NSCs, 
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the evidence suggests that it is an on-going process. The 
nuclei of late-type disk galaxies have been shown to be 
populated by stars spanning a wide range of ages, with 
som e members being as young as a few tens of Myrs 
fe.g. lLong, Charles &: Dubusll2002') . Many authors have sug¬ 
gested that this is indicative of episodic star formation that 
has occurred over an extended period of time. Proponents of 
the GC infall scenario, on the other hand, have argued that 
the young ages of some of the stars can be accounted for if 
star clusters spanning a range of ages were accreted. In at 
least the Milky Way, most young star clusters are of rela¬ 
tively low-m ass in present-day st ar forming regions. On the 
other hand, iNguven et al.l ll2014l l found several super star 
clusters at the centre of Henize 2-10, a blue compact dwarf 
galaxy with an on-going starburst at its centre. The authors 
argue that many of these young massive (~ 10® M©) clus¬ 
ters have very short dynamical friction time-scales, suggest¬ 
ing that this may be a rare snapshot of nuclear star cluster 
formation (via the infall of young massive clusters) around 
a pre-existing super-massive black hole. Regardless, the ev¬ 
idence suggests that GC infall cannot be the whole story in 
at least some NSCs, and that in at least these cases some 
in situ star fo rmation has occu r red. F or example, this was 
illustrated by iHartmann et al.l (1201 il l for the nearby disc 
galaxies NGC 4244 and M33 using both observations and 
dynamical models of the NSCs at their centres. The authors 
argue that gas dissipation is required to account for > 50% 
of these NSCs’ masses. 

In this paper, we calculate and compare the mean- 
square velocities of NSCs and their host galaxies, which we 
use as proxies for energy per unit mass in these systems. 
Traditional observables, including the radial surface bright¬ 
ness and second-order velocity moment profiles, vary on rel¬ 
atively short time-scales due to episodes of matter infall 
and/or star formation. Meanwhile, the total mass, energy 
and angular momentum remain approximately conserved. 
Total energy and angular momentum are decided by the 
various processes that go into forming NSCs and their host 
galaxies, whereas traditional observables are also affected by 
any subsequent or secondary evolution that occurs in energy- 
and angular momentum-space post-formation. Thus, the to¬ 
tal energy and angular momentum more directly probe the 
formation of galaxies and their nuclear star clusters, by of¬ 
fering access to more fundamental properties of the nuclear 
cluster-galaxy system than traditional observables. How¬ 
ever, we caution that our analysis relies on the same observ¬ 
able quantities to calculate mean-square velocities and that 
we make the additional assumption of virial equilibrium, 
which introduces some further uncertainty. Regardless, as 
we will illustrate, the comparison is nonetheless instructive. 
In particular, there is a priori no reason to expect that the 
formation and subsequent time evolution of NSCs should 
be coupled to their host galaxies in energy- and angular 
momentum-space. If NSCs and their hosts form indepen- 
dentljQ, then they should remain independent in energy- 
and angular momentum-space at the present-day. Our re- 


^ For example, while we might naively expect that at least some 
correlation should be present in energy-space in most formation 
scenarios, violent or rapid episodes of mass growth in NSCs could 
erase it. 


suits suggest that NSCs are coupled to their hosts in energy- 
space. We subsequently explore the physical processes that 
could produce such a correlation and, based on this, propose 
a means of also using total angular momentum to further 
constrain NSC formation. In Section[2l we describe the req¬ 
uisite theoretical background. We present the data we use 
in this paper to calculate the root-mean-square velocities in 
Section [3l both for early- and late-type galaxies. Our results 
are presented in Section|31 In Section[3 we discuss the signif¬ 
icance of our results for different NSC formation scenarios, 
and conclude in Section |6l 


2 BACKGROUND 


In this section, we describe the necessary theoretical back¬ 
ground. The purpose here is to motivate our choice for cal¬ 
culating mean-square velocities for both the central NSCs 
and their host galaxies. As we will show, the mean-square 
velocity is a proxy for the energy per unit mass of the sys¬ 
tem, assuming virial equilibrium. The evolution in energy- 
and angular momentum-space is what determines the bal¬ 
ance between the total system mass, half-mass radius and 
root-mean-square velocity (see Equations |4] and [5]), via the 
virial theorem. Thus, the total energy and angular mo¬ 
mentum, which are initially determined by the various for¬ 
mation processes and are subsequently conserved in time, 
are fundamental to determining the observed parameters of 
a self-gravitating system, which are not time-independent. 
Naively, there is no reason to expect that the formation and 
subsequent time evolution of NSCs in energy- and angu¬ 
lar momentum-space should be coupled to their host galax¬ 
ies. If NSCs and their hosts form independently, then the 
evolution in energy- and angular momentum-space should 
proceed independently, first on a crossing time-scale (virial- 
ization) and then on a relaxation time-scale (two-body re¬ 
laxation in energy-space and resonant relaxation in angular 
momentum-space; although the two-body relaxation time- 
scale tends to exceed a Hubble time in both NSCs and, es¬ 
pecially, their host galaxies) (e.g. lMerrittl2013f l. Intriguingly, 
however, we will show that there is a strong (nearly) linear 
correlation, and in fact a near equality, between the energy 
per unit mass (via the mean-square velocity) of NSCs and 
that of their host galaxies. 

We assume virial equilibrium for all NSCs and host 
galaxies. The scalar virial relation yields 2E = -2K = W, 
where E is the total system energy, K is the total kinetic 
energy and W is the total potential energy. The kinetic and 
potential energies are, respectively: 

K = ( 1 ) 

and 


W = - 


GM^ 

rg 


( 2 ) 


where Vrms and rg are the root-mean-square stellar speed 
and the gravitational radius, respectively, and M is the to¬ 
tal stellar mass. It follows from the virial theorem that the 
mean-square velocity is given by: 


_ _ 0.4GM 

M Th 


(3) 
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where the last ( approximate) equality holds for many simple 
stellar systems llSpitzerj|l969l ). and rn is the median radius 
(or half-mass radius) which is related to the gravitational 
radius via rh ~ 0.4rg. In general, the gravitational radius 
is related to the system half-mass radius via a multiplica¬ 
tive constant that depends on the gravitational potential, or 
rg = CTrrh. For example, f or a spheric ally symmetric stellar 
system, we have rg = 2rh djaffdllQ^ ). Throughout this pa¬ 
per, we assume rg = 2rh (i.e. Pr = 2 in Equation [S]) for all 
NSCs as well as for the bodies of all early-type galaxies, and 
rg = rh (i.e. = 1) for the b odies of all late-type galaxies 

(e.g. iBinnev fc Tremaindl 19871) (see below). 

We calculate the mean-square velocities for NSCs and 
their host galaxies using Equations [T] and E quation [2] For 
the NSCs, we use (iBinnev fc Tremaindll98^ : 


^rms,NSC 


GMnsc 

2rNSC 


and for the host galaxies we use: 


^ GMcal 

^r/gr Gal 


(4) 

(5) 


where Mnsc and Mcai are, respectively, the total stellar 
(or baryonic) mass of the NSC and host galaxy, fg is the 
baryonic mass fraction inside roai, and rxsc and roai are 
the effective or projected half-light radii rh of the NSC and 
host galaxy, respectively, which we use as proxies for the 3- 
D half-mass radii (see below). These masses and radii are 
derived from the observations as described in Sections [33] 
and We further assume a constant dark matter mass 
fraction (but see SectionUjfor further discussion of this issue 
and its effects on our results) inside the effective radius rh 
for all gal axies, corresponding to a baryonic mass fract ion 
fg = 0.87 (ICappellari et al.ll2013al f9: iToloba et al.ll2014al f9). 

In practice, kinematical data (i.e. the two terms in 
Equation [6] below, which correspond to the second-order 
velocity moments projected along the line of sight) are 
difficult to obtain. Consequently, very few such mea- 
surements have been provided in the literature (e.g. 
BusajellOj3iOngo&Feol|^ ]j92l : IChae. Bernard! fc KravstovI 
20141 : iToloba et al.ll2014a [bl). This is the reason we calcu¬ 
late root-mean-square velocities indirectly from the virial 
theorem using the total system mass and the half-light or 
effective radius, as opposed to using the directly measured 
second-order velocity moments. 

How does the root-mean-square velocity calculated from 
the virial theorem relate to what is actually observed? Sim¬ 
ilarly to rg, the root-mean-square velocity Vrms is not mea¬ 
sured directly, but rather indirectly via the observed second- 
order velocity moment v^-|- a^, where v and a are the ob¬ 
served mean stellar velocity (i.e. the mean velocity of or¬ 
dered bulk motion) and velocity dispersion, respectively. To 
obtain v^ms from the observed velocity moments projected 
along the ling-o f-sight, we would need to includ e multiplica¬ 
tive factors, or (iBusarello. Longo fc Eeolilll992) : 

2 2.2 

rrma = CtvV + ( 6 ) 


The constants and are determined by the gravi¬ 
tational potential and the angle of inclination relative to 
the line-of-sight. For example, — 1 and = 3 in the 
case of an isothermal s phere, which is spherically symmetric 
llEmsellem et al.ll2007l ). For a non-spherical potential, how¬ 


ever, the angle of inclination must be accounted for when 
measuring the parameters Ov and Oo-. 


3 THE DATA 

In this section, we describe the data used to compile our sam¬ 
ples of both early- and late-type galaxies, beginning with the 
former. We summarize in Appendix|3|all NSC galaxy masses 
and root-mean-square velocities calculated from these data. 
This is done in Tables \M IM and IA3l for each of the late- 
type, Virgo Cluster and Coma Cluster samples, respectively. 


3.1 Early-type galaxies 


For our analysis of early-type galaxies, we use data from 
the Virgo and Coma Cluster Surveys. Beginning with the 
former, we use data for 47 nucleated early-type galax¬ 
ies observed during the Advan ced Camera for Surveys 
Virgo Cluster Survey (ACSVCS. ICote et al.ll20(^. We re¬ 
j ect fi ve galaxies from the original sample of Cote et aP 
(l2006l) in which the apparent NSCs are significantly off¬ 
set from the galaxy’ s phot ocentre and therefore, as dis¬ 
cussed bv ICote et al.l (l2006l ). may well be globular clusters 
that only appear t o reside close to the nucleus due to a 
chance projection llLeigh. Boker fc Kniggel |2012|). We also 
discard four galaxies with extended nuclei ( Graharnll2012l : 
IScott fc Grahaij 20131) . which a re more accurately described 
as nuclear disks ( Balcellsll2007^ . from our sample. 

In order to obtain estimates for the stellar masses, we 
use the apparent z-band magnit udes, (g-z) colors , and half- 
light radii for b oth NSCs (frorn Cote et al]l2006l ) and host 
spheroids (from iFerrarese et al.l 20061 ). IrTorder to convert 
to absolute magnitudes and physi cal radii, we obtain dis- 
tances to indiv i dual g alaxies from iTonrv et al.l (12001[) and 
iBlakeslee et al.l 1 I 2 OO 2 I) wherever po s sible, and for galaxies 
not in the catalogu e of iTon rv et al ] 1 I 2 OOII ) we adopt a dis¬ 
tance of 16.52 Mpc llCote et al.ll200d ). We also need to mul¬ 
tiply the respective z-band luminosities by an appropriate 
mass-to-light ratio. We us e the empir i cally calibrated mass- 
to-light ratios provided by [Bell et al.l (l2003l) . accounting for 
the (g-z) color of NSC and spheroid, respectively. We note 
here that, given the morphological types of the sample galax¬ 
ies (E, SO, dE, dSO, and dE,N), their stellar spheroids can 
be expected to be virialized, and to have little or no current 
star formation activity. This justifies use of a single (colour- 
dependent) mass-to-light ratio for each galaxy spheroid in 
order to derive its stellar mass. Approximate error bars for 
the NSC masses were calculated us ing the 0.041 mag uncer¬ 
tainty quoted bv ICote et al.l ll2006t) . 

We also use data for 53 nucleated low-mass early-type 


galaxies observed durin g the Coma Advanced Camera for 
Surv eys Cluster Survey JCarter et al.l I2OO8I : Iden Brok et al.l 
I2OI4I) . We include all galaxies in the sample designated as 
possible members or better. However, we also check our de¬ 
rived scaling relations without including galaxies listed as 
possible members, since here membership is not certain. Our 
results turn out to be insensitive to the adopted inclusion 
criterion. To calculate the NSC and host galaxy masses, we 
use the corresponding F814W magnitudes. Unlike with the 
Virgo Cluster sample, we do not have colour information 
to calculate empirically-calibrated mass-to-light corrections 
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for the stellar population. Thus, we assume a mass-to-light 
ratio of 2 for all NSCs and host galaxies to correct for the 
unseen component of the stellar mass distribution. This is a 
reasonable assumption for at least the Coma Cluster, since 
with on ly a few exceptions, t he colour distribution is very 
narrow dden Brok et al.ll20l4 l. 


3.2 Late-type galaxies 

For our analysis of late-type galaxies, we compiled a 
sub-sample of 69 galax ies compiled from the sample of 
ICeorgiev fc BokCT (l2014h . for which there are measurements 
of their PEtrosian radii in the 10th SDSS data releas^El- 
These are nearby (< 40 Mpc), mostly low-inclination spirals. 
We calculate the total galaxy luminosity and its effective 
radius using their SDSS Petrosian meas u remen ts following 
the prescription given in iGraham et al.l (l2005ll (their Sec¬ 
tions 3.2 and 3.3, and equations 5 and 6 ). To obtain masses 
for these nuclear star clusters, we use the luminosities mea¬ 
sured from the flux within the best-fitti ng King model of a 
given concentration index (as provided in lCeorgiev fc Bokeil 
ll2014ll l. multiplie d by a suitable mass- t o-ligh t correction ob¬ 
tained using the iBruzual fc Chariot! (l2003l ) simple stellar 
population models for solar metallicity and a lKroupr] ll 200 ll l 
initial mass function. We use the available colour informa¬ 
tion given by the various combinations of the most reliable 
and well-calculated filters, namely F45QW, FSSSW, FGOGW 
and FSldW. All av ailable colours are m a tched to the respec¬ 
tive colours of the IBruzual fc Chariot! ll2003tl models, and 
then a median mass-to-light ratio is calculated for all colours 
(for more details see Georgiev et al. 2015, in preparation). 
For the galaxy luminosities and sizes, we use the z-band 
Petrosian magnitudes and Petrosian effective radii given in 
the tenth data release of the SDSS database, which corre¬ 
spond approximately to the half-light radii of the galaxies. 
The luminosities are converted to masses using the g- and z - 
band Petrosian magnitudes along with the I Bell et al.l (l2003l l 
empirically-calibrated mass-to-light ratio colour-corrections. 


4 RESULTS 

The top panels in Figures [T] [5] and [3] show the relation be¬ 
tween nuclear cluster mass and the mass of its host galaxy 
for all early- and late-type galaxies. The dashed lines have 
slopes of unity and are shifted along the y-axis to approxi¬ 
mately coincide with the distribution of data points in each 
panel. Lines of best-fit are shown for all samples by the solid 
lines, found by performing weighted least-square fits to the 
data. All fit parameters are summarized in Table [T] In all 
samples, the intrinsic dispersion in the data exceeds the un¬ 
certainties on the individual data points. This is partly be¬ 
cause, in some galaxies, the calculated uncertainties are too 
small to be representative of the true underlying uncertain¬ 
ties in the measurements (see Tables lAll IA2I and IA3I in 
Appendix EJ. This is due to the fact that we do not al¬ 
ways have errors for both the masses and radii used in our 
calculations of the root-mean-square velocities. To account 
for the artificially small error bars, uncertainties on the fit 


^ http://cas.sdss.org/drlO/en/home.aspx 


parameters are found by adding (in quadrature) an addi¬ 
tional constant term Uint to the uncertainties on both the 
NSC and host galaxy masses, and forcing the resulting re¬ 
duced chi-square of the fit to be unity. Hence, the term Uint 
is effectively a free parameter in our fits, and is set by the 
condition Xred = 1- This increases the uncertainties on the fit 
parameters to more realistically represent the data. We also 
divide all galaxy masses by 10 ® M©, which is approximately 
equal to the sample means, in order to minimize the uncer¬ 
tainties on the y-intercepts of our lines of best-fit. For the 
Virgo Cluster sample we find a slope of 1.74 ± 0.28 (Virgo), 
whereas for the Coma Cluster and late-type samples we find 
clearly sub-linear slopes of 0.68 ± 0.10 (Coma) and 0.55 
± 0.08 (late-types). The y-intercepts are 6.47 ± 0.07 (late- 
types), 6.17 ± 0.13 (Virgo) and 7.01 ± 0.06 (Coma). For 
comparison, we also perform bootstrapped maximum like¬ 
lihood non-symmetric error-weighted fits (without adding 
an intrisinc dispersion term to the uncertainties on the in¬ 
dividual data points). In all but the late-type sample, the 
derived fit parameters are consistent with our previous esti¬ 
mates to within the uncertainties^ The slopes for the late- 
type, Virgo and Coma samples are, respectively, 0.554 qq|, 
1.921q 53 and 0.7940 25 . The corresponding instrinsic scat¬ 
ter estimates are, respectively, 0 . 21 l 4 ooo 5 , 0.2034Q'gg5 and 
0.19340 005 . T he intrin s ic sca tter and its uncertainty are 
found following I Jeffrey^ (1 19461') and based on Bayesian prob¬ 
ability theory, such that the intrinsic dispersion is invariant 
to rescalings of the problem. 

As seen in Figure[T] the fit is poor for the late-type sam¬ 
ple. This is due to the presence of ~ 10 significant outliers at 
the low galaxy mass end of the distribution, many of which 
have NSCs that are nearly as massive as their host galax¬ 
ies. Hence, these outliers are in part an artifact of having 
only Included the stellar mass In our estimates for the total 
galaxy masses, and not the gas mass. In late-type galax¬ 
ies, this effect can be significant. If the gas mass in these 
outliers is included in the estimates for their total masses, 
these data points shift to the right in Figure [T] and it be¬ 
comes apparent that even these NSCs are only < 1% the 
mass of their host galaxy. Additionally, previous studies have 
shown that, in at least early-type galaxies, the dark matter 
mass fraction decreases with increasing galaxy mass or Ser- 


sic index fe.g.lCappellari et al. 


example, iToloba et al.l (l2014bll recently illustrated that an 


l2006l:lForbes et al.ll2008ll . For 


anti-correlation exists between fpM (within the effective ra¬ 
dius) and total galaxy luminosity for a sample of 39 dwarf 
early-type galaxies in the Virgo Cluster. This effect is not 
included in our fits, and should also contribute to increasing 
the slopes of our lines of best-fit (i.e. closer to unity in the 
Coma Cluster and late-type samples). 

The bottom panels in Figures [T] [5] and [3] show the 
relation between nuclear cluster root-mean-square velocity 
and the root-mean-square velocity of its host galaxy, cal¬ 
culated using Equations U and [SI respectively. As before. 


® The maximum likelihood method used here is the same as in 
Georgiev et al. (2015, in preparation), where it is described in 
more detail. The fit reported here for the late-type sample is 
taken from Georgiev et al. (2015, in preparation), which uses a 
much larger sample of 247 late-type galaxies. This accounts for 
the smaller uncertainties on the fit parameters, compared to the 
Virgo and Coma Cluster samples. 
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Table 1. Parameters for all lines of best-fit for the late-type, Virgo and Coma samples. Each table entry is given in the form (a; /3), 
where a and /3 are the slope and y-intercept, respectively. 


Parameter 

log Mnsc = olog (Mgai/10^) d- /3 

log Vrms,NSC = olog (vrms,gal/40) d- /3 

Late-type 

Virgo 

Coma 

0.55 ± 0.08; 6.47 ± 0.07 

1.74 ± 0.28; 6.17 ± 0.13 

0.68 ± 0.10; 7.01 ± 0.06 

0.85 ± 0.13; 0.16 ± 0.25 

1.34 ± 0.17; -0.78 ± 0.33 

0.58 ± 0.11; 0.61 ± 0.16 


lines of best-fit (solid lines) are obtained from a weighted 
least-squares fit to the data. The uncertainties on the fit pa¬ 
rameters are once again found by accounting for the intrinsic 
dispersion and forcing a reduced chi-square of unity. In the 
Virgo Cluster and late-type samples, the slopes are consis¬ 
tent with being linear to within two standard deviations, or 
1.34 ± 0.17 (Virgo) and 0.85 ± 0.13 (late-type), respectively. 
In the Coma Cluster sample, however, the slope is sub-linear, 
or 0.58 ± 0.11. With these slopes, we find y-intercepts of - 
0.78 ± 0.33 (Virgo), 0.61 ± 0.16 (Coma) and 0.16 ± 0.25 
(late-types). Intriguingly, in both the late-type and Virgo 
Cluster samples, the slopes shift closer to unity relative to 
the corresponding slopes in Figures[T]and [2] Again, for com¬ 
parison, we perform bootstrapped ma:x;imum likelihood non- 
symmetric error-weighted fits, and find that all of these fits 
are consistent with those presented above to within one stan¬ 
dard deviation. For the late-type, Virgo and Coma Cluster 
samples, we find slopes of, respectively, O.QfijlQ ^g, 2.1315 ’q 9 
and l.OSiQ gg. The corresponding instrinsic scatter estimates 
are, respectively, 0.178to:ooL 0 . 182 l ^;225 and 0.182t2;^D^. 
Note that these intrinsic scatter values are comparable to, 
albeit slightly smaller than, those for the corresponding fits 
for the Mgai-MNSC relations discussed above. 

The bottom panels in Figures[T] [^and [3]suggest that, 
to within an order of magnitude, NSCs have roughly the 
same root-mean-square velocity, and hence energy per unit 
mass, as their host galaxies. This is the case despite sub- 
linear slopes for the Coma Cluster sample. For example, the 
fit for the Coma Cluster sample predicts NSCs at low galaxy 
masses that have higher root-mean-square velocities than 
their hosts. However, the distribution of galaxy masses, and 
root-mean-square velocities, is narrow for the Coma Cluster 
sample, so we do not actually observe any NSCs with Vrms 
values much higher than their hosts (i.e. by more than order 
of magnitude). 


5 DISCUSSION 

In this section, we discuss the significance of our results for 
NSC formation. 


5.1 Energy per unit mass 

First, we re-iterate a connection between the root-mean- 
square velocity at the half-light radius and the energy per 
unit mass of the system. From Equations[T]and[21 energy per 
unit mass scales as E/M oc v^^jg oc M/arCg. Therefore, if a 
given NSC has about the same Vrms as its host galaxy (within 
rh), then both systems have approximately the same energy 
per unit mass. As mentioned earlier, this result is highly un¬ 
likely if NSCs form independently of their host galaxies. Our 




Figure 1. The top panel shows the relation between the stellar 
mass of NSCs and that of their host galaxies for our sample of 
late-type galaxies. The bottom panel shows the relation between 
the stellar root-mean-square velocity for NSCs (Equation |4]| and 
the root-mean-square velocity (Equations [SJl of the host galaxy. 
The solid lines show the corresponding lines of best fit given in 
Table [T] The dashed line in the bottom panel shows the one-to- 
one line. 


results are consistent with this general picture, since the fit 
parameters and instrinsic scatter we obtain from our statis¬ 
tical analysis suggest that, in all of the late- and early-type 
samples considered here, the relation between NSC root- 
mean-square velocity and that of its host is as statistically 
significant as the underlying relation between NSC mass and 
host galaxy mass. However, given the uncertainties inherent 
to our analysis, this issue should be re-visited when better 
data become available. 

Why do NSCs have roughly (to within an order of mag¬ 
nitude) the same energy per unit mass as their host galaxy? 
Said another way, our results illustrate that the average star 
in a nuclear star cluster has roughly the same kinetic energy 
(and hence total energy, by the virial theorem) as the aver¬ 
age non-cluster star in the host galaxy. One interpretation 
of this result is that NSCs were effectively formed from stars 
taken directly from the host galaxy, while conserving energy. 
Hence, in this picture, galaxies form first, and then their 
NSCs subsequently form from some small piece of their host 
galaxy bodies. From our perspective, this scenario seems to 
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Figure 2. Same as Figure^ but for the Virgo Cluster sample. 




log ^^. 0.1 (in km/s) 


Figure 3. Same as Figure[T]but for the Coma Cluster sample. 


be the simplest physically-motivated picture that is consis¬ 
tent with our results. However, it is by no means the only 
Interpretation. For example, it is conceivable that NSCs and 
their host galaxies form with different root-mean-square ve¬ 
locities, and then somehow evolve over time toward a com¬ 
mon root-mean-square velocity. However, we are unaware 
of any physical mechanism that might contribute to such a 
trend. 


5.2 Implications for NSC formation and evolution 


Our result can be understood within the framework of both 
the GC infall and in situ star formation models for NSC 
formation. For example, consider the GC infall model for 
NSC formation in a host galaxy with constant circular ve¬ 
locity at all radii (i.e. both inside and outside the effective 
radius rh), an idealized assumption we will return to below. 
The key point is that an infalling GC (on roughly a circular 
orbit) that has just reached the nucleus due to dynamical 
friction will have approximately the same velocity relative 
to the centre of the galaxy as it did several kiloparsecs out. 
Hence, when stars orbiting within the GC have reached the 
galaxy centre, their velocities (relative to the galaxy centre 
of mass) remain similar to the circular velocity at or be¬ 
yond rg. It follows that, when the inspiraling GC reaches 
the galaxy’s centre and becomes part of the central nuclear 
cluster, the energy coming from the bulk motion of the in¬ 
falling GC is converted to internal thermal energy within 
the nuclear cluster itself. That is, for a stationary observer 
at the centre of the galaxy, the centre of mass of the orbit¬ 
ing cluster transitions from having a net velocity roughly 
equal to the circular velocity to having zero net velocity, 
whereas the stars within the cluster always orbit the galaxy 
centre of mass at the circular velocity. Thus, for galaxies 
with a constant circular velocity inside the effective radius, 
this predicts a linear relationship between the energy per 
unit mass in the central NSC and that of its host galaxy 
(or, equivalently, a linear relation between the root-mean- 
square velocity of the central NSC and that of its host, at 
all radii inside rg). 

Next, consider the in situ star formation scenario for 
NSC formation, again in a host galaxy with constant circu¬ 
lar velocity at all radii. Ignoring radiation losses, we would 
naively expect a similar picture as outlined for the GC infall 
model. That is, the gas should reach the central NSC with 
the same velocity as it had several kiloparsecs out, creat¬ 
ing a disk (assuming a preferred plane of accretion) of star¬ 
forming gas with a circular velocity that roughly matches 
that of its host (ignoring any central SMBH). We caution, 
however, that the picture is potentially more complicated 
within the framework of the in situ star formation model 
for NSC formation. This is because, when infalling gas fila¬ 
ments collide with the forming nucleus, shocks can develop 
that will radiate energy away and lower the energy per unit 
mass of the central NSC relative to its host. Thus, at least in 
the ideal case of a collisionless star-forming gas, we expect 
a central NSC with roughly the same energy per unit mass 
as its host in galaxies with a constant Vc at all radii, just as 
in the GC infall model. 


The assumption of a constant circular velocity as a 
function of galactocentric radius is, to first-order, consis¬ 
tent with the observations for most galaxies in our sam¬ 
ples (ignoring the innermost regions near the galactic cen¬ 
tre). First, consider the early-type galaxies. Observations 
have now revealed that early-type galaxies are composed 
of both a spherical bulge and an underlying disk (e.g. 


Emselleni_e£_^ 2007: Graham fc WorlevI l2008l : Scott et al.l 
20141 : Laurikainen et al.ll2011^ . The ratio of bulge-to-disk lu- 
minosity in typical early-type galaxies is Lb /Lp ~ 1/3 (e.g. 
ICraharn fc Wori^l200^ : iLaurikainen et al.ll201ll l. Thus, for 
our early-type samples, the potentials can be approximated 
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by an isotherm al density profile for the bu lge, combined with 
a Mestel disk dBinnev &: Tremaine! 1 198m . Subsequently, in¬ 
tegrating over the surface brightness profile yields a mass- 
dependence of the form M(r) oc r. The key features of 
the se potentials is that they have circular velocities given 
by (iBinnev fc Tremainelll987| '): 


GM{r) 


(7) 


where r is the distance from the centre of the galaxy, and 
M(r) is the mass distribution within r. The surface bright¬ 
ness is proportional to 1/r in both potentials, hence their 
sum also scales as 1/r. Subsequently, integrating over the 
surface brightness profile yields a mass-dependence of the 
form M(r) oc r. From Equation [71 it follows that the circu¬ 
lar velocity is independent of r in the combined bulge-disk 
potential, provided the outer extent of the bulge extends 
beyond rh. Beyond the limiting radius of the bulge, it con¬ 
tributes a constant to the total potential, and Vc declines 
weakly with increasing r. For our sample of late-type galax¬ 
ies, the potential can, to first-order, be approximated by a 
simple Meste l disk, for which the circu lar velocity is inde¬ 
pendent of r (iBinnev fc Tremainelll98^ 1. Thus, the circular 
velocity is independent of r in this potential as well. There¬ 
fore, in both early- and late-type galaxies, the assumption 
of a constant circular velocity Vc within the effective radius 
rh is a reasonable first-order assumption. 

We caution that the potentials of real galaxies are 
considerably more complicated than the assumptions made 
above suggest. In particular, the ratio M(r)/r drops consid¬ 
erably at small galactocentric radii, roughly where the NSC 
potential begins to dominate over that of its host. Previous 
studies have used numerical simulations to model the prop¬ 
erties of NSCs forme d from GC infall (e.g. lAntonlnl et al.l 
20121: lAntoi^l2013ll . focusing on the final few parsecs. 


Antoninil ( 2011 1 found that, after many mergers of GCs 


with the central NSC, the NSC mass-radius relation steepens 
from Rnsc ~ M^lc to Rnsc Mnsc- Thus, the energy per 
unit mass in the NSC ends up being independent of the num¬ 
ber of accreted GCs, and rea ches a roughly constant value. 
Using more detailed models, ICnedin, Ostriker fc Tremaine! 
ll2014h found a similarly weak dependence of the form Rnsc 
~ Mnsc- Therefore, our results suggest a similarly weak de¬ 
pendence of the galaxy half-mass radius on the total galaxy 
stellar mass, in order to reproduce the observed relation be¬ 
tween the NSC root-mean-square velocity and that of its 
host. Of equal importance, the circular velocity profiles of 
real galaxies often are not constant at large galactocen- 
tric radii, but inste ad show considerable fluctuations (e.g. 
iToloba et al.|[^14alf9l . Thus, our assumption of a constant 
circular velocity at all radii is not strictly correct, and is 
likely the source of much of the scatter seen in Figures [T1 [2] 
and O 

What else might a linear relation in energy per unit 
mass be telling us about NSC formation and/or evolution? 
In particular, we ask: How might a linear relation in energy 
per unit mass between NSCs and their host galaxies, once in 
place, be affected by other physical processes characteristic 
of galaxy formation, such as collisions and mergers occur¬ 
ring in galaxy clusters? Within the framework of the above 
picture for NSC formation, which should approximately pro¬ 
duce NSCs with the same energy per unit mass as their 


host, one possible interpretation of this result is that any 
collisional events affecting the energy per unit mass of the 
host must have similarly affected its NSC. In other words, 
significant energy gains or losses cannot have occurred for 
either the NSC or its host galaxy post-formation. This is be¬ 
cause additional energy deposited within NSCs and/or their 
host galaxies due to direct collisions could contribute to in¬ 
creasing the half-light radius of the system, and hence to 
decreasing the energy per unit mass (provided the particles 
Involved in the collision are themselves collisionless, which 
is a decent assumption for stars). This is because the ki¬ 
netic energies of the colliding galaxies are positive and, via 
energy conservation, must be added to the total energy per 
unit mass of the collision/merger product. This contributes 
to a reduction in the total energy per unit mass of the sys¬ 
tem (which is a negative quantity for bound systems), which 
might be observed as an increase in the half-light or effective 
radius. Thus, if NSC formation involved more energetic colli¬ 
sions than occurred in the host galaxy, this could contribute 
to a lower root-mean-square velocity in the NSC relative 
to its host galaxy. Conversely, if galaxy formation involved 
more energetic collisions than the formation of the central 
NSC, this could contribute to a lower Vrms in the host rel¬ 
ative to its central NSC. If the efficiency of either of these 
scenarios scales with total galaxy mass, then this could con¬ 
tribute to a sub-linear or even super-linear dependence of 
NSC root-mean-square velocity on host galaxy Vrms- In gen¬ 
eral, collisions that affect the NSC more (or less) than its 
host contribute to a higher (or lower) energy per unit mass 
in the host relative to the central NSC, and should increase 
the scatter in the bottom panels of Figures [T1 [2] and [S] The 
intrinsic scatter in the observed relations is at least consis¬ 
tent with being due in part to different merger histories for 
the galaxies in our samples. 

We caution that the presence of a central SMBH can 
also strongly affect the evolution and observed structural 
parameters of NSCs. For example, the presence of a cen¬ 
tral NSC might be a required pre -cursor to the forma¬ 
tion of an SMBH via core collapse llMiller fc Davie j|2012l : 
iGnedin. Ostriker fc Tremainell2014l l. Subsequently, the cen¬ 
tral SMBH can ac t to disrupt infalling clusters before they 
reach the nucleus llAntonini et al.l[2012 : lAntoninill2013ll . Al¬ 
ternatively, a central SMBH or SMBH-SMBH binary could 
contribute to a higher energy per unit mass in NSGs rela¬ 
tive to their hosts, by ejecting stars with a large (positive) 
energy per unit mass. More work is needed to properly quan¬ 
tify the impact of the presence a central SMBH on the re¬ 
sults presented in this paper. Finally, we have also not con¬ 
sidered any internal evolution that might occur in NSCs, 
affecting their structural parameters, nor have we consid¬ 
ered energy exchange betw een an existing NSC and its host 
galaxy (e.g. lMerritdl2013h . In general, these omissions are 
reasonable since the relevant time-scales (e.g. two-body re- 
laxation) ten d to exceed a Hubble time in most NSCs (e.g. 
lMerrittl[2^1 . 

5.3 Future work 

We have argued 1) that NSCs form after their host galaxy 
bodies, and 2) that, in principle, both formation scenarios 
for NSC formation discussed above can produce NSCs with 
the same energy per unit mass as their host galaxies. In 


© 2011 RAS, MNRAS OOO.IlIfTTl 






































8 Leigh N. W. C., Georgiev I. Y., Boker T., Knigge C., den Brok M. 


an attempt to identify the most plausible mechanism for 
NSC formation, we now turn our attention to angular mo¬ 
mentum. The general predictions discussed below that arise 
from consideration of angular momentum could be testable 
in the future with higher resolution surveys. 

What do we expect for the angular momentum content 
of NSCs relative to their host galaxies? We will argue that, 
unlike energy per unit mass, the answer could depend on the 
NSC formation mechanism. In particular, GC orbits form a 
halo within the host galaxy, with their orbital planes aligned 
at random angles relative to each other. Thus, if only a sin¬ 
gle GC inspiral is responsible for forming the NSC, then we 
would expect significant rotation in the remnant NSC, with 
the axis of rotation aligned perpendicular to the original or¬ 
bital plane of the inspiraling GC. This rotation can be can¬ 
celed by the cumulative effects of many GCs inspiralling into 
the nucleus, which should create a central NSC with roughly 
the same energy per unit mass as its host, but nearly zero 
net angular momentum (within the framework of the pic¬ 
ture described above) Q This is the case for both early- and 
late-type galaxies. The story changes within the framework 
of the in situ star formation model for NSC growth. This is 
because the cumulative effects of galaxy mergers have not 
depleted the total angular momentum content of late-type 
galaxies. In fate-types, there is a preferred plane of accretion 
for gas, but not for any accreted GCs. Therefore, any in situ 
star formation occurring in a NSC forming via gas accretion 
should generally occur in a disk, with its dominant angu¬ 
lar momentum vector aligned with that of the host. NSCs 
formed from many episodes of GC infall, however, should 
have roughly zero net angular momentum. Thus, we predict 
that NSCs in late-type galaxies formed primarily via gas ac¬ 
cretion and in situ star formation should contain significant 
angular momentum, which is directly observable as rotation. 
Naively, this also predicts a correlation between the observed 
inclination of the host galaxy disk relative to the observer 
line of sight and the (isophotal) ellipticity of the NSC, with 
more edge-on disk galaxies having more flattened NSCs. In 
early-type galaxies, on the other hand, there is no preferred 
plane of accretion for the gas. Thus, we expect a central NSC 
with roughly the same energy per unit mass as its host, just 
as in the GC infall model, but with little to no net angular 
momentum or rotation. 

The implication that NSCs are formed and evolve to¬ 
gether with their host galaxies suggests that the origins of 
NSCs can be distinguished observationally (again, ignoring 
any internal evolution occuring within the NSCs) when bet¬ 
ter resolution becomes available. This can be done by de¬ 
composing NSCs into bulge and disk components, and com¬ 
paring the bulge and disk masses, as well as the angle of 
inclination between the plane of the NSC disk and that of 
its host galaxy. This is especially true in spiral galaxies, since 
they have a milder merger history than ellipticals. Consider 
the quantity P = cos(9Mdisk/(Mbuige+Mdisk), where Mdisk is 

In fact, repeated episodes of GC infall may still yield non- 
negligible rotation in the NSC, provided the orbits characteris¬ 
tic of the bulk of the GC population in the host galaxy exhibit 
a clear rotation signature, with v/a > 0.5. This has been ob¬ 
served in a numb er of early-type galaxies fe.g. lBeaslev et ali gOOgt 
Pota^^l l2013^■ and even the bulge of the Milky Wav fe.g. lHar^ 
lOOh^^OloTi^ateh . 


the mass of the disk component, Mbuige is the mass of the 
bugle component and 9 is the angle of inclination between 
the NSC disk and that of its host galaxy. If /? < 0, then 
the NSC was formed mainly through GC infall. If, on the 
other hand, /3 > 0 then a non-negligible fraction of the NSC 
formed through in situ star formation, with (3 = 1 corre¬ 
sponding to the case of 0% GC infall and 100% in situ star 
formation. 


6 SUMMARY 

We estimate and compare the root-mean-square velocities 
for a large sample of nuclear star clusters and their host 
galaxies (both early- and late-type). These are used as prox¬ 
ies for energy per unit mass, and it is demonstrated that 
NSCs have roughly the same energy per unit mass as their 
host galaxies, to within an order of magnitude. The origin 
of this interesting relation is discussed. We interpret this as 
evidence that NSCs do not form independently of their host 
galaxies, but rather that their formation and subsequent 
evolution are coupled. We discuss how our results can po¬ 
tentially be used to offer a clear and observationally testable 
prediction to distinguish between the different nuclear star 
cluster formation scenarios, and even quantify their relative 
contributions. 


APPENDIX A: DATA 

In this appendix, we present the data used to generate Fig¬ 
ures [T1 |2]and [3l All calculated NSC and galaxy masses and 
root-mean-square velocities are summarized in Tables \M\ 
IA2I and IA3I for each of the late-type, Virgo Cluster and 
Coma Cluster samples, respectively. 
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Table Al. Properties of all late-type galaxies. Column 1 gives the object ID. Columns 2 and 3 list the total stellar mass of the galaxy 
and NSC in units of 10^ Mq and 10® Mq, respectively. Columns 4 and 5 provide the root-mean-square velocities for the galaxy and 
NSC, respectively, in units of km s“^. The galaxy and NSC mass uncertainties are taken directly from the literature. The uncertainties 
for the root-mean-square velocities are calculated using the provided uncertainties on both the galaxy or NSC masses and radii, wherever 
available. 


ID Mgai Mnsc ^rms,gal Vj-mSjNSC 

(10® M0) (10 ® Mq) (km s-l) (km s"!) 


8.76l5!;g®o®gg 7.963 ± 1.71 

0.127'^;gip% 1.239 ± 0.266 


IC0769 

M108 

MCG-01-03-085 

NGC0428 

NGG0600 

NGG0853 

NGG0864 

NGG1042 

NGG2500 

NGG2541 

NGG2552 

NGG2805 

NGG3041 

NGG3259 

NGG3274 

NGG3319 

NGG3338 

NGG3346 

NGG3359 

NGG3423 

NGG3445 

NGG3455 

NGG3666 

NGG3756 

NGG3913 

NGG3949 

NGG4030 

NGG4041 

NGG4062 

NGG4096 

NGG4204 

NGG4208 

NGG4237 

NGG4393 

NGG4395 

NGG4411b 

NGG4498 

NGG4517 

NGG4522 

NGG4525 

NGG4534 

NGG4567 

NGG4571 

NGG4592 

NGG4618 

NGG4625 

NGG4631 

NGG4635 

NGG4651 

NGG4701 

NGG4900 

NGG4904 

NGG5112 

NGG5204 

NGG5300 

NGG5334 

NGG5585 

NGG5668 

NGG5964 
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0 0190.003 
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U.OZ‘±_Q QQ 4 

0.2335 ± 0.0153 
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15.03 ± 0.5514 
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zi. /U-o 216 
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3-4520j30|83 

4.066 ± 0.6458 

U.U10_Q QQ4 

3.648 ± 0.4912 

0.592005®, 

8.186 ± 1.758 


78.12 ± 0.7954 

69.95 ± 

15.02 

65.68 ± 3.675 

32.65 ± 

7.009 

42.96 ± 0.8615 

42.34 ± 

7.954 

63.53 ± 3.214 

80.54 ± 

28.87 

68.75 ± 0.5458 

65.52 ± 

15.3 

69.64 ± 0.6369 

52.51 ± 

12.26 

105.1 ± 4.838 

249.3 ± 

6.466 

35.41 ± 11.38 

85.62 ± 

64.64 

39.13 ± 1.374 

19.6 ± 

5.41 

41.26 ± 1.527 

36.04 ± 

7.738 

24.06 ± 1.237 

87.37 ± 

19.93 

77.24 ± 1.515 

103.7 ± 

23.67 

105.9 ± 4.083 

56.52 ± 

26.09 

136.6 ± 0.8513 

72.56 ± 

15.27 

47.78 ± 2.064 

24.41 ± 

2.293 

45.37 ± 0.688 

28.37 ± 

4.505 

138.3 ± 3.337 

198.5 ± 

20.37 

72.56 ± 0.4986 

77.74 ± 

18.16 

101.1 ± 0.2946 

54.31 ± 

6.618 

70.02 ± 0.9634 

41.84 ± 

3.823 

60.51 ± 0.2824 

75.91 ± 

17.32 

91.52 ± 0.5246 

89.05 ± 

19.12 

125.6 ± 0.9413 

104.1 ± 

62.37 

120.7 ± 0.5633 

158.4 ± 

86.31 

58.8 ± 0.6037 

10.67 ± 

2.45 

142.9 ± 0.1635 

101.9 ± 

23.29 

283.4 ± 10.26 

268.1 ± 

13.64 

220.6 ± 0.4428 

93.96 ± 

55.21 

131.7 ± 0.2014 

52.71 ± 

11.32 

108.5 ± 1.624 

32.75 ± 

7.032 

22.89 ± 0.7458 

21.39 ± 

4.996 

141. ± 0.4547 

39.25 ± 

8.427 

174.5 ± 0.5642 

68.42 ± 

10.68 

3.289 ± 136. 

19.05 ± 

3.852 

11.53 ± 1.301 

67.77 ± 

13.33 

38.07 ± 0.7635 

46.21 ± 

10.54 

25.77 ± 3.636 

34.06 ± 

7.313 

27.99 ± 4.157 

28.84 ± 

6.191 

77.42 ± 1.023 

17.49 ± 

3.756 

66.77 ± 0.3455 

64.58 ± 

12.77 

46.78 ± 0.5342 

20.89 ± 

3.91 

106.8 ± 2.657 

238.6 ± 

50.21 

91.27 ± 1.382 

83.46 ± 

2.165 

56.28 ± 0.2334 

39.64 ± 

8.51 

62.28 ± 0.3859 

15.66 ± 

8.069 

58.72 ± 0.2616 

20.69 ± 

0.699 

21.62 ± 2.386 

17.57 ± 

3.696 

48.05 ± 0.719 

47.38 ± 

17.61 

251.2 ± 0.721 

172.1 ± 39620000, 

114.9 ± 0.665 

218.4 ± 49.82 

106.4 ± 5.072 

113.6 ± 3.698 

95.14 ± 0.857 

44.3 ± 10.35 

50.58 ± 1.697 

21.43 ± 1.306 

34.09 ± 0.4149 

35.43 ± 2.32 

68.9 ± 2.784 

39.61 ± 1.453 

93.79 ± 0.77 

39.5 ± 3.276 

47.78 ± 0.964 

34.58 ± 5.492 

71.93 ± 5.838 

54.91 ± 8.721 

20.7 ± 2.778 

57.18 ± 7.698 

114.2 ± 6.918 

33.59 ± 7.212 




10 Leigh N. W. C., Georgiev I. Y., Baker T., Knigge C., den Brok M. 

Table A2. Properties of all Virgo Cluster galaxies. The columns are the same as in Table ixn 


Cappellari 


ID Mgai 

(10® Mq) 


VCC1720 

20.21 

± 5.31 

VCC1883 

16.66 

- ± 3.92 

VCC1242 

15.27 

' ± 3.17 

VCC784 

16.79 

i ± 3.14 

VCC828 

13.69 

i ± 3.00 

VCC1250 

3.68 

± 2.04 

VCC1125 

8.06 

; ± 3.4 

VCC1283 

9.08 

± 1.79 

VCC1261 

4.87 

± 1.69 

VCC698 

9.52 

± 2.56 

VCC1422 

3.82 

± 1.35 

VCC2048 

2.94 

± 1.08 

VCC1871 

2.26 

± 0.58 

VCC1910 

2.1 

± 0.81 

VCC856 

2.22 

± 0.85 

VCC140 

2.33 

± 0.71 

VCC1355 

1.82 

± 0.35 

VCC1087 

3.29 

± 1.07 

VCC1861 

2.88 

± 1.00 

VCC543 

2.19 

± 0.65 

VCC1431 

2.2 

± 0.74 

VCC1528 

1.63 

± 0.48 

VCC1695 

1.69 

± 0.78 

VCC437 

2.8 

± 1.25 

VCC2019 

1.02 

± 0.72 

VCC33 

0.43 

± 0.25 

VCC200 

1.34 

± 0.49 

VCC1488 

0.41 

± 0.38 

VCC1895 

0.76 

± 0.39 

VCC1545 

1.41 

± 0.42 

VCC1192 

1.85 

± 0.62 

VCC1075 

1.11 

± 0.47 

VCC1627 

1.03 

± 0.32 

VCC1440 

1.19 

± 0.44 

VCC230 

0.69 

± 0.29 

VCC2050 

0.29 

± 0.26 

VCC751 

1.41 

± 0.51 

VCC1828 

1.26 

± 0.47 

VCC538 

0.63 

± 0.26 

VCC1407 

1.24 

± 0.42 

VCC1886 

0.51 

± 0.1 

VCC1199 

0.58 

± 0.16 

VCC1539 

0.52 

± 0.11 

VCC1185 

0.89 

± 0.51 

VCC1826 

0.61 

± 0.22 

VCC1489 

0.32 

± 0.06 

VCC1661 

0.43 

± 0.08 


M. et al. 2006, MNRAS, 366, 1126 


MnSC 
(106 Mo) 


89.24 ± 3.371 

25.12 ± 0.949 
23.14 ± 0.8741 
63.61 ± 2.403 

104.4 ± 3.942 

21.48 ± 0.8115 
3.814 ± 0.144 
8.99 ± 0.3396 

14.25 ± 0.5383 
11.33 ± 0.4278 
7.342 ± 0.2773 
2.105 ± 0.0795 

24.18 ± 0.9134 

10.1 ± 0.3814 

19.66 ± 0.7427 
0.7763 ± 0.0293 

2.38 ± 0.0899 
8.816 ± 0.333 
6.08 ± 0.2296 
1.074 ± 0.0406 

10.49 ± 0.3963 
0.7213 ± 0.0272 

1.06 ± 0.0400 

6.18 ± 0.2334 
5.671 ± 0.2142 
0.7211 ± 0.0272 
0.392 ± 0.0148 
0.1285 ± 0.0049 
0.2053 ± 0.0078 
1.146 ± 0.0433 
19.78 ± 0.747 
2.178 ± 0.0823 
38.08 ± 1.438 
14. ± 0.5287 

5.445 ± 0.2057 
0.6627 ± 0.0250 
2.203 ± 0.0832 
1.567 ± 0.0592 
2.382 ± 0.0900 
4.252 ± 0.1606 
0.9599 ± 0.0363 
14.29 ± 0.5396 
3.209 ± 0.1212 

2.6 ± 0.0982 
7.802 ± 0.2947 
0.5612 ± 0.0212 
5.142 ± 0.1942 


Vrms,gal 
(km s“^) 


142.9 ± 37.55 

143.4 ± 33.73 

168.1 ± 34.91 

210.1 ± 39.29 

176.3 ± 38.64 

96.18 ± 53.32 

149.3 ± 62.98 
120. ± 23.66 

86.41 ± 29.99 
139. ± 37.39 

73.24 ± 25.88 

84.73 ± 31.12 

106.3 ± 27.29 
73.46 ± 28.34 
64.05 ± 24.53 
88.08 ± 26.84 

42.9 ± 8.25 

70.6 ± 22.96 

69.81 ± 24.24 
60.79 ± 18.04 

82.82 ± 27.86 
71.36 ± 21.01 
43.91 ± 20.27 

57.66 ± 25.74 

40.64 ± 28.69 
38.06 ± 22.13 
56.15 ± 20.53 
34.72 ± 32.18 

47.42 ± 24.34 

63.66 ± 18.96 

89.74 ± 30.08 
43.87 ± 18.58 
90.63 ± 28.16 
69.07 ± 25.54 

48.38 ± 20.33 
27.61 ± 24.76 

57.9 ± 20.94 
48.56 ± 18.11 

63.25 ± 26.1 
59.28 ± 20.08 
34.27 ± 6.72 

73.43 ± 20.26 
26.78 ± 5.665 
36.03 ± 20.65 
52.35 ± 18.88 
29.84 ± 5.596 

18.38 ± 3.42 


Vrms.NSC 
(km s“^) 


158.5 ± 14.36 

167.7 ± 49.31 

133.2 ± 27.12 

111.1 ± 6.401 
102. ± 5.162 

145.9 ± 39.67 
41.31 ± 5.066 

67.49 ± 9.271 

103.1 ± 20.42 

86.12 ± 15.06 
75.05 ± 15.28 

42.7 ± 9.777 
82.14 ± 6.162 
84.47 ± 15.88 

58.74 ± 3.485 
30.78 ± 9.863 
41.01 ± 7.712 

93.64 ± 24.53 
37.04 ± 2.589 

12.13 ± 0.631 
34.77 ± 1.678 

32.8 ± 12.82 

34.44 ± 10.13 
44.71 ± 4.132 
72.46 ± 17.7 
24.6 ± 5.461 

16.64 ± 3.13 

11.75 ± 3.319 

16.2 ± 5.435 
28.84 ± 5.563 

66.25 ± 4.577 
38.73 ± 7.103 
72.04 ± 3.736 
81.93 ± 10.7 
66.56 ± 14.34 
16.42 ± 1.849 

41.12 ± 8.369 

27.17 ± 3.491 

46.17 ± 10.91 
29.98 ± 2.003 
25.07 ± 4.384 
78.03 ± 9.158 
18.03 ± 0.831 
37.36 ± 5.418 
91.54 ± 25.86 

16.12 ± 2.038 
41.03 ± 3.83 


R., Naab T., Oosterloo T., Sarzi M., Scott N., Serra P., 


Cappellari M., Scott N., Alatalo K., Blitz L., Bois M., Weijmans A.-M., Young L. M. 2013b, MNRAS, 432, 1862 

Bournaud F., Bureau M., Crocker A. F., Davies R. L., Capuzzo-Dolcetta R., Miocchi P. 2008, MNRAS, 388, 69 


Davis T. A., de Zeeuw P. T., Due P.-A., Emsellem E., Carter D., et al. 2008, ApJS, 176, 424 


Khochfar S., Krajnovic D., Kuntschner H., McDermind 
R. M., Morganti R., Naab T., Oosterloo T., Sarzi M., 
Serra P., Weijmans A.-M., Young L. M. 2013a, MNRAS, 
432, 1709 

Cappellari M., McDermind R. M., Alatalo K., Blitz L., 
Bois M., Bournaud F., Bureau M., Crocker A. F., Davies 
R. L., Davis T. A., de Zeeuw P. T., Due P.-A., Emsellem 
E., Khochfar S., Krajnovic D., Kuntschner H., Morganti 


Chae K.-H., Bernard! M., Kravstov A. V. 2014, MNRAS, 
in press (arXiv: 1305.5471) 

Cote P., Blakeslee J. P., Ferrarese L., Jordan A., Mei S., 
Merritt D., Milosavljevic M., Peng E. W., Tonry J. L., 
West M. J. 2004, ApJS, 153, 223 
Cote P., Piatek S., Ferrarese L., Jordan A., Merritt D., 
Peng E. W., Hasegan M., Blakeslee J. P., Mei S., West M. 
J., Milosavljevic M., Tonry J. L. 2006, ApJS, 165, 57 
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Table A3. Properties of all Coma Cluster galaxies. The columns are the same as in Table ED Uncertainties on the NSC masses are 
provided when available. 


ID 

Mgal 

(10® Mo) 

MnSC 
(10® Mq) 

VrmSjgal 

(km s“^) 

Vrms,NSC 
(km s“^) 

LEDA126789 

2.339 ± 0.021 

12.28 ± - 

79.96 ± 0.19 

30.76 ± 

1.10 

SDSSJ125950.18-275445.4 

2.133 ± 0.019 

8.492 ± - 

57.47 ± 0.20 

37.2 ± 

1.69 

SDSSJ130007.12-275551.4 

1.476 ± 0.013 

25.65 ± - 

73.16 ± 0.52 

48.86 ± 

1.48 

SDSSJ125926.45-275124.7 

1.346 ± 0.012 

2.812 ± 0.259 

58.06 ± 0.29 

20.71 ± 

2.35 

SDSSJ130026.16-280032.0 

1.12 ± 0.010 

8.492 ± - 

47.36 ± 0.29 

40.03 ± 

1.40 

SDSSJ125914.43-280217.3 

0.162 

± 

0.006 

3.707 ± - 

16.2 ± 0.43 

21.1 ± 

1.42 

SDSSJ125953.93-275813.7 

1.476 

± 

0.013 

9.312 ± - 

55.6 ± 0.17 

32.9 ± 

1.25 

SDSSJ125636.78-271247.8 

1.618 

± 

0.028 

3.707 ± - 

40.25 ± 0.60 

28.97 ± 

3.05 

SDSS J130000.97-275929.5 

1.021 

± 

0.009 

21.33 ± - 

47.33 ± 0.30 

54.4 ± 

1.40 

SDSS J125844.58-274458.2 

1.476 

± 

0.013 

17.74 ± - 

44.28 ± 0.10 

39.71 ± 

1.15 

SDSSJ130042.86-280313.8 

0.849 

± 

0.007 

2.339 ± - 

34.21 ± 0.19 

23.51 ± 

2.84 

COMAil25949.960p275433 

0.588 

± 

0.005 

10.21 ± - 

43.51 ± 0.30 

37.28 ± 

1.42 

SDSS J130032.61-280331.4 

0.706 

± 

0.012 

11.2 ± - 

30.2 ± 0.52 

39.41 ± 

1.27 

SDSSJ130036.58-275552.2 

0.588 

± 

0.005 

4.887 ± - 

36.76 ± 

0.31 

31.79 ± 

2.14 

COMAil25713.240p272437 

0.588 

± 

0.005 

7.064 ± - 

37.44 ± 

0.18 

30.16 ± 

1.26 

SDSSJ125942.36-280158.5 

0.644 

± 

0.011 

5.358 ± - 

31.93 ± 

0.42 

20.43 ± 

0.96 

SDSS J130027.57-280323.9 

0.536 

± 

0.009 

1.618 ± - 

25.91 ± 

0.40 

17.55 ± 

1.86 

SDSSJ130004.03-280030.7 

0.489 

± 

0.004 

7.064 ± - 

36.71 ± 

0.43 

33.8 ± 

1.53 

SDSS J125902.43-280021.3 

0.407 

± 

0.004 

5.875 ± - 

36.47 ± 

0.27 

29.62 ± 

1.64 

SDSSJ130018.70-275512.6 

0.371 

± 

0.003 

4.065 ± - 

31.06 ± 

0.22 

29.42 ± 

1.45 

SDSSJ130044.10-280215.4 

0.536 

± 

0.009 

5.875 ± - 

25.56 ± 

0.47 

28.55 ± 

1.10 

SDSS J125943.53-275620.6 

0.371 

± 

0.010 

9.312 ± - 

26.72 ± 

0.54 

42.67 ± 

1.55 

COMAil25828.358p271315 

0.489 

± 

0.009 

3.381 ± - 

27.16 ± 

0.33 

24.12 ± 

1.35 

SDSSJ130042.51-280325.4 

0.338 

± 

0.006 

7.745 ± - 

32.48 ± 

0.53 

27.71 ± 

0.89 

SDSS J130037.30-275441.0 

0.338 

± 

0.009 

8.492 ± - 

27.2 ± 

0.68 

38.22 ± 

1.53 

SDSSJ125955.93-275748.6 

0.281 

± 

0.002 

2.565 ± - 

25.61 ± 

0.30 

16.91 ± 

1.31 

SDSSJ130032.96-275406.6 

0.536 

± 

0.009 

11.2 ± - 

21.68 ± 

0.42 

39.41 ± 

1.27 

SDSS J130003.18-275648.3 

0.148 

± 

0.003 

1.476 ± 0.413 

18.94 ± 

0.28 

15.75 ± 

2.95 

SDSSJ125927.22-275257.0 

0.195 

± 

0.003 

5.875 ± - 

21.92 ± 

0.30 

18.71 ± 

0.62 

SDSS J125951.46-275935.4 

0.213 

± 

0.004 

3.707 ± - 

23.99 ± 

0.33 

24.04 ± 

1.22 

SDSS J125930.83-275810.2 

0.234 

± 

0.014 

7.745 ± - 

23.96 ± 

1.30 

35.52 ± 

1.35 

SDSSJ125945.55-280313.4 

0.162 

± 

0.007 

3.381 ± - 

17.93 ± 

0.56 

17.78 ± 

1.31 

COMAil3035.990p275505.4 

0.162 

± 

0.003 

2.565 ± - 

18.36 ± 

0.31 

22.19 ± 

2.18 

SDSS J130024.85-27.99085 

0.162 

± 

0.004 

6.442 ± - 

20.25 ± 

0.46 

31.02 ± 

1.29 

SDSS J125853.08-274741.8 

0.234 

± 

0.008 

3.083 ± - 

22.85 ± 

0.61 

22.94 ± 

1.46 

SDSS J125708.35-272923.9 

0.148 

± 

0.005 

5.875 ± - 

22.63 ± 

0.75 

31.06 ± 

1.42 

SDSS J125959.08-275841.4 

0.281 

± 

0.007 

11.2 ± - 

18.53 ± 

0.62 

43.89 ± 

1.40 

SDSSJ130030.94-280312.8 

0.085 

± 

0.004 

3.707 ± - 

14.98 ± 

0.65 

24.3 ± 

1.44 

SDSSJ130030.94-280312.8 

0.093 

± 

0.005 

3.707 ± - 

14.86 ± 

0.76 

23.13 ± 

1.24 

COMAil25937.351p28210.6 

0.085 

± 

0.003 

1.618 ± 0.149 

17.2 ± 

0.59 

14.02 ± 

1.66 

SDSS J130034.32-275817.6 

0.093 

± 

0.002 

1.476 ± 0.136 

16.91 ± 0.30 

15. ± 1.81 

SDSSJ125934.39-275942.9 

0.148 

± 

0.004 

1.945 ± - 

19.09 ± 0.46 

16.17 ± 1.31 

SDSS J 130039.32-275748.4 

0.093 

± 

0.003 

2.133 ± - 

14.9 ± 0.43 

16.88 ± 1.15 

COMAil3039.554p275350.0 

0.123 

± 

0.003 

2.565 ± - 

17.57 ± 0.49 

22.92 ± 2.18 

SDSSJ125951.81-275726.3 

0.135 

± 

0.012 

5.875 ± - 

15.34 ± 1.44 

31.18 ± 1.68 

COMAil25925.477p28211.0 

0.071 

± 

0.006 

1.476 ± 0.136 

13.8 ± 0.97 

11.72 ± 2.03 

SDSSJ125700.89-273155.1 

0.102 

± 

0.011 

2.339 ± - 

11.54 ± 0.88 

21.01 ± 1.47 

SDSSJ130030.94-280312.8 

0.085 

± 

0.004 

3.707 ± - 

14.98 ± 0.65 

24.3 ± 1.44 

SDSSJ130030.94-280312.8 

0.093 

± 

0.005 

3.707 ± - 

14.86 ± 0.76 

23.13 ± 1.24 

SDSS J130029.81-280401.0 

0.085 

± 

0.002 

1.021 ± 0.094 

14.54 ± 0.31 

12.18 ± 1.81 

SDSSJ130004.04-275342.7 

0.281 

± 

0.015 

5.875 ± - 

14.38 ± 0.72 

30.94 ± 1.41 

SDSS J125952.18-275946.3 

0.078 

± 

0.002 

1.228 ± 0.113 

15.05 ± 0.40 

15.15 ± 2.24 

SDSSJ125856.78-274644.5 

0.049 

± 

0.001 

0.234 ± 0.088 

13.76 ± 0.31 

3.34 ± 2.08 
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